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r'-J ■ Abstract 

><■ 

We study at finite temperature the energy-momentum tensor T liv {x) 
of (i) a scalar field in arbitrary dimension, and (ii) a spinor field in 1+1 
dimensions, interacting with a static background electromagnetic field. 
T^ v separates into an UV divergent part T^ a representing the virtual 
sea, and an UV finite part Tf^ sma describing the thermal plasma of 
the matter field. T^f" remains uniform in the presence of a uniform 
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electric field E, while T^ sma becomes a periodic function with pe- 
riod Ax = 2-kT/E in the direction parallel to E. A related periodicity 
is found for a uniform static magnetic field if one spatial direction 
perpendicular to the magnetic field is compactified to a circle. 
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1 Introduction 



A finite-temperature (T > 0) or thermal quantum field can be visualized as 
a sea of virtual particles occupying space together with a thermal gas of real 
field excitations. The virtual particle sea is independent of the temperature 
T. It can, however, be deformed by coupling the field to static background 
structures of various kinds. This is generally known as the static vacuum 
Casimir effect. Somewhat less widely known is that, for boundaries and 
other static backgrounds, the thermal gas is "mechanically" distorted along 
with the sea. 

Abelian (and non-Abelian) gauge theories present Casimir problems with 
particular features arising from the underlying gauge invariance. Indeed, 
the restrictions on the class of allowed gauge transformations imposed by 
T > are found to have remarkable consequences for the spatial energy 
distribution of a charged thermal matter field coupled to a static background 
electromagnetic field (E,B). The local distortion by (E,B) of the virtual 
sea and thermal plasma (as we now refer to the thermal gas consisting of 
both particles and antiparticles) can be revealed by a local analysis in terms 
of the thermal stress energy momentum tensor. 

The problem of charged quantum fields coupled to a uniform electromag- 
netic background field is an old one, going back to famous papers of Euler 
and Heisenberg M and Schwinger 0. There has been done a great deal of 
subsequent research, much of it reviewed in refs. 0-0. However, global as- 
pects of the problem have received most of the attention, while local aspects 
seem to have been neglected. 

In the present paper we try to gain general insight into the local response 
of gauged thermal matter fields to a background electromagnetic field. Using 
the Matsubara formalism (see e.g. refs. [0, || ||), we work on a hypercylinder 
of circumference (3 = 1/T in the euclidean time direction, choosing space to 
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be flat and infinite. Our results reveal some surprising features: We find 

— * 

for both scalar and spinor fields that a uniform background electric field E 
causes the thermal plasma to become non-uniform - in fact periodic in its 
spatial distribution along the direction of the E, with period Ax = 2irT/\E\. 
The sea nevertheless remains spatially uniform. 

The organization of this paper is as follows: We begin by considering ther- 
mal scalar fields, for which the discussion can easily be carried out in arbitrary 
dimension. By coupling the scalar field to an arbitrary static gauge potential 
A^(x) we show that the characteristic effects arising from the minimal cou- 
pling are common to all dimensions. Special attention is drawn to periodicity 
features related to gauge invariance and to the topology S 1 xR d of space-time. 
We then compute for the specific gauge potential = (Exi + const, 0) the 
thermal stress tensor Tf£\x). The periodicity of the thermal plasma along 
x\ with period Axi = 27rT/E is thereby made explicit. 

In section 3 we turn to the Schwinger model at T > with the linear 
background potential referred to above. The discussion is now more involved, 
but the calculations reveal features similar to those encountered in the scalar 
case. We conclude in section 4 by commenting on a claim in the litera- 
ture concerning the factorization of the thermal heat kernel, and drawing an 
analogy between our findings and the Quantum Hall Effect. 

2 Thermal scalar field 

Scalar electrodynamics is useful as a theoretical laboratory for studying gauge 
theory phenomena in arbitrary space-time dimension. We first consider the 
general problem of a thermal scalar field coupled to an arbitrary static i) 
background potential V(x) and ii) gauge potential A^x). We then specialize 
to the potential A^ = (Exi + const, 0) for a uniform background electric field 
E = (E,0, ...,0) and compute explicitly the thermal stress tensor of <p. The 
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case of a uniform magnetic field is also discussed briefly. 



2.1 Scalar field in a static background Schrodinger po- 



To set the stage we briefly review the case of a scalar quantum field interacting 
with a static background potential V(x) in d- dimensional free space R d . We 
wish to study the thermodynamical properties and vacuum Casimir energy of 
this system. To this end it will be convenient to work in the imaginary time or 
Matsubara formalism. Euclidean space-time is then a hyper- cylinder S 1 x R d . 
Correspondingly we impose periodic boundary conditions in euclidean time 
on the scalar field 4>(xq, x) 



where (3 = 1/T and T is the temperature. The spectral operator for the 
theory in question is [— <9q — A + V(x)} with A the Laplacian in d dimensions. 
The vacuum and thermodynamical properties of the system can be computed 
from the bilocal heat kernel 



where Af and (f) k (x ,x) are the eigenvalues and respective eigenf unctions of 
the spectral operator, 



tential 



4>(x ,x) = 0(x o + 



(2.1) 



h^(t;x,y)=J2e^MxW k (y) 



k 



[-dl-/\ + V{x)]<p k {x)=\l<p k {x) . 



With subject to the boundary condition ( [2.1|) we have 



(j) k (x) — > (fi mn (x) 




'mn 




(2.2) 



4 Throughout this section we use the euclidean notation x = {x^} = (xq,x). 
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where the spatial modes (p n (x) and associated spectum {a>^} are determined 
by the spatial mode equation 

[-A + V(x))(p n (x) = Lul<f n (x) . 

The Euclidean thermal Green function f\ then has the spectral representation 

<n*my)>0 = L rr^rz 7^1 

m,n [\ ft I ~r "-riJ 

J0 P m n 

We may perform the Matsubara sum by using the Jacobi identity [p"0| 

£ e -6(m-a)» = £ e -^e- jW (2.3) 



, 6 

m=— oo Z=— 00 



with the result 



< U x )Uv) >P = / dth(t;x,y) T=0 





00 22 

x E e-^e^ (x °- yo) (2.4) 



i=— 00 

where 



h{t;x,y) T=0 = 44 x 

n 

is the T = bilocal heat kernel of the operator [— <9q — A + V(a?)]. Hence, for 
a static background potential the T > Green function separates |TT| into 
two distinct and well defined parts: the virtual sea part (1 = contribution) 
which is independent of T and coincides with the T = Green function, and 
the thermal gas part (I 7^ contribution) which exhibits the full temperature 
dependence and vanishes exponentially as T — > 0: 



< (j){x)(j){y) >=< <f){x)<f){y) > sea + < <t>{x)<t>{y) > 



gas 



3 We denote operators by a " hat" . 
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Now suppose that we have calculated from < 4>(x)4>{y) > sea the T = 
vacuum stress tensor 

TZ=<f^>=Y,T,r (2.6) 

n 

as a (still to be renormalized) spatial mode sum. Repeating the calculation 
at finite temperature we obtain from eq. fl2.4| ) 



1 _|_ p-P^n 



rpflU _|_ 

sea ^ / y n 



o 0UJ n 

n ^ 



rpiw _|_ a-i/i!/ /o 7\ 

sea ^ gas ■ V 1 / 



At finite temperature the mode sum is thus modified by the familiar Bose- 
Einstein distribution, in agreement with one's expectations. This completes 
our brief review of a usual Casimir problem at T > 0. 

2.2 Scalar field in a static background gauge potential 

Let us now consider a (massive or massless) scalar quantum field (f>(x) coupled 
to a static background Abelian gauge potential A^x). We shall again be 
interested in studying the thermodynamical properties and Casimir energy 
of this system. 

On a space-time hyper- cylinder S 1 x R d , both <fi(xo,x) and A^xq^x) are 
required to satisfy periodic boundary conditions in euclidean time: 

(f){x ,x) = (f)(x + f3,x) ,A^(x ,x) = A IM (x + P,x) (2.8) 

The relevant spectral operator in this case is the gauged Laplacian in d + 1 
dimensions —D 2 , where = <9 M — iA^ couples the quantum scalar field to 
a static background gauge potential A^x). (We have absorbed the electric 
charge into A^.) The thermodynamical and vacuum properties of the system 
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can again be computed from the bilocal heat kernel 

hW{t;x,y) = J2e-^e~ tM2 M^my) 

k 

where A^ are the eigenvalues of —D 2 , M is the mass and <f>(x) is subject to 
the boundary condition ( |2.8| ). Periodicity in x implies 

<Pk{x) -> 4> mn = -^e* 2 7 a ' B Vmn(z) , 

where ip mn (x) now satisfies the associated eigenvalue problem 
-(V -iA(x)) 2 + V m (x) 

with 

V m {x) = [A (x) - ^y-f . (2.10) 

Notice that the m-dependence of the Schrodinger-like background poten- 
tial V m (x) leads to a coupling of spatial position x with the Matsubara fre- 
quencies. Eqs. (|2~9l ) thus represent a different equation for each Matsubara 
frequency ^p, with n labelling the complete set of normalizable solutions 
{(Pmn} of this Schrodinger problem for a given potential V m (x). The situ- 
ation is thus very different from the scalar case discussed previously. It is 
characteristic of gauge theories and has very important consequences as we 
shall see. 

On the cylinder S 1 x R d we can always gauge a static Aq(x) to the interval 
[0, H\ , but not in general to zero, if we respect the periodicity property (|2.8|) 
of A^ and <fi. The only exception is when A = N 2 ^-, in which case we 
may gauge Aq to zero by performing the (allowed) gauge transformation 
A — A + d \ with A = (^f)x . 

Another way of stating this is to observe that in the exceptional case 
Aq = N^j the gauge transformation can be absorbed into the Matsubara 
index m via the transformation m —>■ m + N. For this reason Aq(x) is always 
gauge equivalent to a configuration taking values in the range [0, ^]. In 
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<Pmn(x) = \ mn <fmn (2-9) 



the zero-temperature limit, on the other hand, we may always gauge Aq(x) 
to zero. Indeed, the discrete Matsubara frequencies k Q = become a 
continuous variable k in the range — oo < k < oo, and A (x) may be 
absorbed into a shift in k under the integral J dk . 

Green function 

Following the steps of the previous section, we have this time for the 
thermal Green function 



r dt e ' m2 \ E e 2 -T (x °- yo) h m (t; x, y) (2.11) 
Jo p m 



where 

h m (t-x,y) =Y, e ^ X2mn( Pmn(xW* mn (y) (2.12) 

n 

is the spatial bilocal heat kernel for the m'th spatial background potential 
V m (x) in (|2.10| ) . Notice that because of the coupling of Matsubara frequen- 



cies with the spatial modes the Green function can no longer be straight- 
forwardly separated into "sea" and "gas" contributions as in the potential 
problem discussed previously. In special cases this separation can neverthe- 
less be made explicit, as we shall see. 

Energy-momentum tensor 

The symmetric canonical energy-momentum tensor for the complex scalar 
field (f)(x) coupled to a background is formally given by (in Minkowski 
space-time) 



where 



t= l -[{D^{D^)-M 2 ^}. 
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Using the equation of motion {D^D^ + M^cf) = we have for the divergence 



of Tpy, 



and for the trace, 

f; = - l -{d- l)(Djy(D^) + I(d+ l)M 2 0t0 . (2.13) 

From ( p.!3| ) we see that is traceless only in d = 1 spatial dimension for 
M — 0. We furthermore see that T^q is conserved if F i0 = 0, that is, if 
the electric field vanishes. This still allows for a static magnetic field, which 
makes sense, since a static magnetic field cannot do work on charges. 

From ( |2.4| ) the thermal stress energy tensor Tj£> =< T^ v >@ is now easily 
written down. We have for the separate components 

r°° 9 1 97T 

< J D o 0(.x)[ J D o 0(x)]t> /3 = / dte- tM \{^YY\™-<Z)?h m {t^,y) 

Jo p p m 

(2.14) 

<A^)[D 3 ^)] t > /3 = / dte~ tM -J2^S(D^D]h m (t;x : y)} 

Jo p m x ^y 

(2.15) 

<^(,)[D J #)]t>, = j( dte~ tM2 -{— )£[m-a(£)] 

x limipf^M^y)} (2-16) 

roo 1 

> = / dte- tM -J2h m (t;x,x) (2.17) 
Jo ™ 



< \Mx)\ 2 



where 



Z7T 



is the rescaled temporal component A . 

The above expressions giving the energy-momentum tensor in terms of the 
heat kernel need of course to be properly UV-regularized. We see that these 
results differ essentially from the nongauge results given earlier, since in the 

10 



present case is expressed as a Matsubara sum over nonidentical spatial 
problems, the latter depending on Matsubara label m. Thus an explicit 
separation into a sea and thermal gas (or more properly, thermal plasma) 
contributions is in general not possible. We now consider specific background 
potentials for which the problem is completely soluble, and the separation 
into gas and thermal plasma contributions can be displayed. 

2.2.1 Constant gauge potential 

To begin with we consider a constant background gauge potential 

27T 

A = — a, A = 0. (2.18) 

As already pointed out, on the cylinder S* 1 x R d this gauge potential cannot 
be gauged to zero; however, it is gauge equivalent to a potential with a in the 

— * 

range < a < 1. On the other hand, the spatial component A of a constant 
A^ can always be gauged to zero on this cylinder, so that we may choose 
1=0. 

Following our general notation we have in this case (for infinite volume 
the index n becomes the continuous momentum label k) 

<t>mk{x) = -^—lZ ik '* 

(2tv) 2 

X mk = (j) 2 (m-af + P (2.19) 



and 



V m (x) = (j) 2 (m-af 



h m (t;x,y) = e-^^h^x-y) (2.20) 



where 



h (t;x-y) = J—Jd d ke- tp e rk <^ 



(2tt) c 

1 (3-y) 2 

d e ** 

(47Tt)# 

11 



is the infinite volume, zero temperature heat kernel of the free scalar field. 

The factorization of the heat kernel h m into an m-dependent and and m- 
independent factor now enables one to perform Matsubara sums explicitly 
by making use of the identity 



1 oo 

- -L x - _/ _„,-' 

2b 



e 



-oo 



^ OO 22 

2(-)i n 2 e~^- cos(2vrcm) (2.21) 

" n=l 



obtained from (|2.3| ) by differentiation with respect to b. One finds for the 
Minkowskian [] energy density after some simple algebra 

t-i(/3) rpsea , rpplasma /c\ c\c\\ 

J 00 _ J 00 00 I^Z.ZZJ 

where the temperature-independent part representing the sea is given by the 
UV-divergent integral 

Too a = lfdte- tM \4nt)-^i^ + M*] , (2.23) 

and the temperature dependent part representing the thermal plasma carries 
all the dependence on the gauge potential and is finite: 



rpplasma = / ^ ^ ^ -(^) g 

Jo ~, 



oo 

'"iran) 

2/32 



r n 2 B 2 d-1 n 2 /3 2 

x { i^" + ^r +M }e ~^- (2 - 24) 



The integral in (|2.24f) can be evaluated in terms of the modified Bessel func- 
tion K v (z), but we shall not do so. It is important to observe that the thermal 
part vanishes exponentially as T — > 0, and that it exhibits the expected peri- 
odicity property in the parameter a, in line with our earlier observation that 
a can always be chosen to lie in the interval [0, 1]. 

The remaining diagonal stress tensor elements Ta(f3) with i = l,2,..,d 

are given by Tqq(P) with the signs of the (d — l)/2£ and M 2 terms in the 
6 Note that the transition from Euclidean to Minkowskian space-time involves a minus 
sign:< -|A)0| 2 >p->< \D <t>\ 2 > fi . 
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curly brackets of eqs. ( |2.23| ), (|2.24| ) reversed. All nondiagonal elements of 
vanish. This tensor is obviously conserved. 
Finally observe that for d = 1 and M = we have T^ a = 0, and there is 
no need to perform an UV renormahzation. rurthermore ig = Id 



in this case, in accordance with the trace condition 



2.2.2 Constant electric field 

Next we consider the linear background potential 

A (xi) = Exi + 27ca/(3, A = (2.25) 

corresponding to a constant background electric field £ = (-E, 0, 0) in the 
Xi direction. Here the constant term 2ira/f3 in A has the obvious physical 
significance of determining the location X\ = —2ira/(3E of the zero of this 
linear potential. Global aspects of what now follows have been investigated 
previously for scalar fields as well as for spinor fields (see e.g. refs. ||-|!). 
We focus here on local aspects of this problem. Note that the limit E — > is 
distinctly nontrivial. For that reason we have chosen to analyse separately 
the £^0 problem here and the E = problem in subsection 2.2.1 above. 

Continuing with our general notation (where now n — >• (n, kj_)) we have 
the spatial modes 

= ^(xJ^Trr^V^ (2.26) 

with k± = (k 2 ,...,k d ) and x± = (x 2 , Xa) representing momentum and 
position perpendicular to x\. Inserting <p mn % into the spatial mode equation 
L9|) one obtains 



d 2 2 2 



¥n{Xm) = e n ip n (x m ), (2.27) 



where 



2ty 

x m =xx + — (a-m). (2.28) 



13 



This is just the harmonic oscillator eigenvalue problem in Schrodinger theory 
with orthonormal eigenfunctions 



n\ V n 

e n = 2E(n+-), n = 0,l,2,... (2.29) 

Here H n (z) are Hermite polynomials satisfying y" — 2zy' + 2ny = 0. In eq. 
( p.9| ) the m-dependent backgrounds V m (xi) = E 2 x 2 n are identical harmonic 
oscillator potentials centered at equidistant positions x\ — (m — a)2ir//3E. 
As we shall see this periodic arrangement of identical potentials leads to a 
periodic structure along x\ (with period Axi = 2tt//3E) in T^ u {x) and in 
other local quantum functions. Notice that the m dependence of V 9 mn fc ± (^) 
resides entirely in the argument x m of the harmonic oscillator wave function 
(Pn{xm) - i.e. entirely in the position x\ — (m — a)2n/i3E of the zero of 
V m (xi). One consequence of this fact is that the spectrum of — D 2 given by 

X 2 mnk± = 2E(n + \) + k\ (2.30) 

does not depend on the Matsubara label m. This is in sharp contrast with 
the nongauge scalar theory of section 2.1 with its spectrum ( p. 2D , and with 
the constant problem above with spectrum ( p. 19 ). 



The spatial heat kernel ( 2.12f) constructed from the spatial modes (|2.26|) 

is 



h m (t;x,y) = k E (t;x m ,y m )h (t; \x± - y±\) (2.31) 



where 



oo 

\2 



k E {t;x m ,y m ) = J2 e tK fn{x m )<f* n {y r< 

n=0 



E 



2 

e -i£'(x'i-j/i) 2 coth(2Et) 



27rsinh(2M) 

y^ e -Ex rn y m ta,nh(Et) (2 32) 

and h is the free-space heat kernel for d — 1 dimensions. The mode sum J2 n 
here has been performed with the help of the identity (see e. g. ref. [[RJ, p. 
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194) 



E 

n=0 



n\ \2 



H n (x)H n (y) = (l-z 2 )- 1 / 2 exp 



■[2xy(l - z) - z(x-y) 2 } 



Alternatively one can use known formulae for the propagator in the harmonic 
oscillator problem (see section 3). Note that in the limit E — > the heat 
kernel (|2.32|) smoothly becomes the one-dimensional free-space heat kernel 
as it should 



k E {t 



1 



47T 



e -(a;i-yi) 2 /4t £ 



0. 



The limit — ^ is nonetheless far from being uniform: As E —>■ 
the potentials V m (xi) = E 2 x 2 n — > (2irm/j3) 2 change into m-dependent con- 
stants (like mass terms) independent of x±. The background returns to the 
constant potential = (27ia//3,0) of the preceeding subsection with its 
spectrum A^ fc — (rn — a) 2 (2n/l3) 2 + k 2 . Remarkably, the dependence on m 
so conspicuously absent from the E > spectrum X 2 nnk± reenters the E = 
spectrum. 



The Green function ( 2.11|) and local quantities derived from it possess 
the sea + plasma structure one expects to find. Let us display this for the 
Green function ( [2.11|) in the limit y — > x; 

2Et 



< 



x r >/ 



/■oo „ 

/ dte- tM (Ant) 
Jo 



-d/2 



sinh(2£t) 



x ■ 



-Ex^ tanh(.Ei) 



^2.33) 



where, using the identity (|2.3| ), the Matsubara sum can be evaluated with 
the result 

E 



-i oo 
" m=— oo 



Exi, tanh Et 



Atx tanh Et 

x e -in(2Tra+xif3E) 



E e 

n=— oo 



-n 2 P z E/4 tanh Et 



[2.U) 



Thus we find 



< Mx)\* > f 
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<|0(x)| 2 > 



2 ^plasma 



d-1 
2 . 



E 



dte (47rt) 

An sinh i^t 

oo 

]T 2cosn(27ra + a; 1 /3£)e- n2 ^ /4tanhi? ' 



;2.36) 



n=l 



One easily verifies that < 



>sea coincides with the corresponding T = 



quantity < 



>t=o as it should. This function of course needs UV 



renormalization. All dependence on temperature is in < 



2 ^plasma 



I > 







The 



latter function is finite and it vanishes exponentially as T — > 0. Moreover 
it is periodic in X\ with period Axi = 2ir//3E, reflecting the equidistant 
arrangement of potentials V m (xi) = E 2 x 2 m . 

We now proceed to the straightforward calculation of Tj$). Using eqs. 
( p,14| )- p.l7] ) one easily verifies 

(2.37) 



^oo ( x ) 



rpsea , rpplasma / \ 
1 00 "T J 00 \ X l) 



where 



1 r°° o 



4(d-i). 



E 



4-7T sinh Et 



d- 1 
2t 



+ M 2 



(2.38) 



and 



rpplasma / \ 
-'OO 



o 4tt sinh 

oo 



X E e 

n=\ 



X 



n 2 /3 2 S/4tanh £t 

d-1 



+ M z + n l 



cosn(2na + x\f3E) 
( E(3 \ 2 ' 



(2.39) 



2t \2smhEt) 
The energy densities of the virtual sea and thermal plasma have been ob- 
tained by performing the Matsubara sums in eq. (|2.37|) with the help of 
the identities (|2.3| ), fl2.21| ). For the reader's convenience we give the form in 
which the latter identity is used here: 



rp2 oo 

^ V J -x 2 m Eta,nhEt 

r J TXl— — OO 



E 



Air tanh Et 



E 



x 



2 tanh Et 

oo 

I ^ e -n 2 /3 2 E/4ta,nhEt e -in(2Tra+x 1 l3E) ^ 4Q) 

n=— oo 
o2 p oo 

Y / n 2 e- n2p2E/4tanhEt cosn((3Ex 1 +27Ta)} 



tanh 



n=l 
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Of course, T s q q needs UV renormalization. In the limit E — - > 0, T S Q fl and 
rppiasma ^ove smoothly become the E = functions ( 2.23|) , ( 2.24 ). 

The other diagonal components of T JJ 3 ' can be similarly obtained. For 
brevity we do not write them down. Again T x ± is given by eqs. ( |2.37|) 
- (|2.39| ) with the signs of the (d — l)/2t and M 2 terms reversed. Thus Tqq q = 
_ rpsea £ or an y S p a ^j a } dimension d > 1. All nondiagonal components of T 
vanish. Interestingly T^ a still vanishes for a massless scalar field in d = 1 
spatial dimensions. Indeed the tracelessness of implies T S Q a = T(f a and 
this, combined with T S Q a = —Tff a , leads to T™ a = 0. However, for d > 1 
this tensor does not vanish. 

As expected Tqq° is independent of position: the uniform electric field 
leaves the virtual sea spatially uniform. Physically this seems reasonable. 
Virtual particles do not have the prolonged existence needed to participate 
in e.g. thermal equilibrium: the sea remains temperature- independent. 

Things are quite different for the thermal plasma. The particles of the 
thermal plasma do have prolonged existence, and they do participate in ther- 
mal equilibrium. Also they are nonuniformly redistributed by the background 
electric field E as eq. (|2.36|) shows. The plasma becomes spatially nonuni- 
form - in fact periodic as already described with period Axi = 2-kT/E. As 
temperature T — ► with E fixed, the plasma ceases to exist and spatial 
uniformity is restored. Alternatively, holding T fixed and letting E — > 0, the 
period Axi = 2n/(3E becomes infinite, thereby restoring spatial uniformity 
to the plasma. In fact, the total energy contained in the plasma changes 
discontinuously as we let E tend to zero. Indeed, if we integrate TQQ asma (x) 
over an integral number of periods the result is zero. Hence, the total energy 
of the thermal plasma vanishes in the presence of the uniform electric field. 
However, in the limit E — > the nonzero energy density (|2.24|) characterizes 
the uniform plasma throughout space. 
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2.2.3 Constant magnetic field 

To investigate the effect of a uniform background magnetic field on the virtual 
sea and thermal gas of a scalar field it is of particular interest to consider the 
case d = 3 spatial dimensions. We choose for the static background potential 
= (0, 0, Bxi, 0) leading to the magnetic field B = (0, 0, B). Following our 
general notation we then have for the spatial modes (now n — *> (n, k 2 , k 3 )) 

In 

where k 2 , k 3 are continuous momentum labels in finite space and 

k 2 

Xk 2 =Xi + — . 



The modes <p n {xk 2 ) are the harmonic oscillator wavefunctions ( p. 29 ) satisfy- 
ing 

-df + B 2 x\^ (f) n (x k2 ) = 2B(n + -)0„(x fc2 ) . 
The eigenvalues of — D 2 are now 

and are independent of k 2 - 

It is already apparent that the uniform magnetic field does not introduce 
spatial non-uniformity into either the virtual sea or the thermal gas. Indeed, 
all mode sums involve an integration in k 2 over the infinite interval [—00, 00]. 
Since k 2 is a continuous variable we can perform the shift k 2 — > k 2 — Bx± 
in the integration variable, thereby absorbing the Xi-dependence into the 
integration. Thus T^ v and other local quantum functions will not depend on 

Xi. 

If, however, we compactify the x 2 direction perpendicular to the magnetic 
field to a circle of perimeter L, we are led to a problem paralleling the one 
with constant electric field discussed in the previous subsection. Compact 
x 2 corresponds to discrete momenta k 2 = p(^), with p running over all 
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integers (as in the case of the Matsubara index). The harmonic oscillator 
mode equation above becomes 

-d{ + B 2 x 2 p ip n (x p ) = 2B(n + ^)(p n (x p ) , 

where x p = x\ + p^r- Again we have an infinite set of harmonic-oscillator 
potentials equally spaced at intervals Axi = along the xi-axis. Much as 
in section 2.2.2, local quantities such as T^ v are periodic in x\ with period 
Axi- Clearly one is observing here something akin to the Quantum Hall 
Effect. 

One could of course extend this discussion to include backgrounds with 
both constant electric and magnetic fields. We shall not persue these inter- 
esting matters here, but rather turn to the case of massless fermions in 1+1 
dimensional space-time. 



3 The finite temperature Schwinger model 

Electrodynamics of massless fermions in 1+1 dimensions (the Schwinger 
Model (SM)) is described by the Lagrangian density 

C(x) = --F^{x)F^(x) + $(x) (ift + e4) if>(x). (3.1) 

Our discussion will follow the general lines of the previous section. Aside 
from the usual complications of spin, we shall have to face in the infinite 
volume limit the existence of an infinite number of normalizable zero modes. 



3.1 The case of no electric field 



We choose again Q2.18D for our gauge field configuration. We emphasize 



once more that at finite temperature the Aq field cannot be gauged to zero. 
Hence Aq = const implies observable effects []. The corresponding eigenvalue 
7 The same applies to A\ if we were to compactify space to a circle S . 

19 



equation for the Dirac field may be written in the form ^ 



2-n 

id + —a) + 7 5 <9 1 



(3.2) 



with the solution 



1 



, , - i(2m+l)%x° 

Wm,k{ X ) = — e 3 



1 Pk(Xl) 

\ (PkM 



(3.3) 



where 



(Pk{xi) 



2tt 



A 



and the eigenvalues 



where 



\ m (k) = ±Jal + k2 



, n .x 71 " 27T 



(3.4) 



and 



A = W;rrf A 



respectively. We normalize the eigenfunctions by choosing \A\ = 4^. Note 
that zero modes are absent, unless a = (M + |). 

Three local quantities are of primary interest to us: The heat kernel 
and corresponding zeta-function densities, as well as the external field Green 
function for the construction of the energy-momentum tensor. 

Local heat kernel and zeta-function 

From its definition we have for the matrix valued heat kernel 



h w (t;x,y) a ,p = (x\[e- 1 



{ipf 



alk 



(3.5) 



Our euclidean conventions are 70 = — ci, ,71 = 02 ,7 5 = 03 
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with the corresponding expression for the local zeta function, defined by 

1 



C(s;x,y) 



T(s) 



(3.6) 



Since the eigenvalues come in pairs of positive and negative numbers, the 
"chiral" parts of the 2x2 matrix cancel, leaving us with a diagonal spinor 
structure: 



1 



• (3-7) 



V° 1 



We shall be interested only in the diagonal part, x = y of the heat kernel. 
Taking the trace we have 



trh^(t;x,y)= 1 -^J^e-^) 



(3.8) 



Making use of the Jacobi identity (|2.3|) we may cast tr/i^(t;x,x) into the 
form 



trh^ (t; x, x) = trh(t; x, x) 



1 + 2 ^ (-l) m cos(2vrma)e- m2 « 



m=l 



(3.9) 



where h(t; x, x) is the zero-temperature diagonal heat kernel 

1 r dk 



tih(t; x, x) 



-tk 1 



Indeed, for (5 — > oo, 



-far, 



^TTtJ 27T 



dk° 



(3.10) 



_ (3 J 2tt V47rt 

so that we are lead smoothly from (|3.9|) to ( p,10| ) in the limit /3 — > oo. 

The heat kernel (|3.9|) is interesting in two respects: i) It exhibits the 
remarkable periodicity property in Aq already encountered in section 2. As 
commented in that section, this periodicity property has its origin in the 
fact that Aq(x) is always gauge equivalent to a configuration taking values in 
[0, ii) This exact result provides a counterexample to a claim |l2j made 
in the literature concerning the factorization of the temperature dependence 
in the thermal heat kernel of a gauged spinor field. 
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Energy density 

The energy density is given by 

Too(x) = tr [t l0 D (x)G(x, y; A)] x=y (3.11) 

where G(x, y; A) a p is the external field Green function having spectral rep- 
resentation 

G(x,y;AU = £ J ^ ^ (3.12) 



with the sum extending over the positive and negative eigenvalues 
From Q3.12j ) and fl3.3| ) we obtain for the Green function (recall that the eigen- 



values occur in pairs (A, —A), which results in the off-diagonal structure) 

P m J 27T |A w (fcJ| \ a m -ik 

\ V CLm+ik 



(3.13) 



Performing the differentiation in (|3.11| ) and combining terms, we obtain 



Notice that this sum diverges and hence requires regularization. For the 
special case of a vanishing gauge-potential (a = 0) the result fl3.14| ) has a 
familiar interpretation. Making use of the identity 

^ 1 7T TTX 

/ -, 77: = — tanh — 

^ (2m + l) 2 + x 2 2x 2 



7Tt=— OO 



we may write Q3.14Q for = in the form 



W = |E/£-/ h- 2 / ^W^t ■ (3-15) 

where we have set \k\ = u(k). The last term evidently corresponds to the 
usal Fermi-distribution for a gas of massless fermions, while the second term 
can be interpreted as the negative energy contribution of the Dirac sea. As 
for the first term, one can identify it with 5 2 (0), which is independent of /?; 
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hence this term is eliminated by temperature-independent renormalization, 
as should be the case. 

Effective action and thermal energy 

The effective (euclidean) action is given by 

Seff = ~\ In det {%$ +e4) 2 = ±C'(0) (3.16) 
where ((s) is the zeta-function defined for the spectrum Q3.4T) , 

«"=S/i(5«k (3 ' 17) 

and where the "prime" on indicates differentiation with respect to the 
argument. With the identification S e ff = —InZ, with Z the external field 
partition function, the average thermodynamic energy is thus given by 

In order to allow for a comparison with previous results, we first perform the 
differentiations and then take the limit s — > in the resulting expression, 
which then takes the form 

( 9 . A r dk 1 d . 



It is now important to observe that our parametrization of the temporal part 
of the gauge field in terms of A = ^-a correctly exhibits the /^-dependence 
of the eigenvalue spectrum. Indeed, as we already remarked above, A Q may 
always be mapped into the interval [0, by an allowed gauge transforma- 
tion, and correspondingly the parameter a is taken as lying in the interval 

[0,1]. With this observation we have 

(9 2 

u \ 2 / j,\ _ __ 2 

d/3 m{) ~ ~f m ' 



so that finally 

c/=--c(o) = -r/- al 

in agreement with ( |3.14| ). These considerations show that it is appropriate 



to define the divergent sum (p. 14 ) in terms of the zeta-function via ( |3.17| ). 
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3.2 A constant Electric Field 



We now turn to the case of a constant electric field £ = —E/e, with the 



choice ( |2.25D for the gauge potential. The eigenvalue equation ( |3.2|) is now 
replaced by 



27T 

[ido + Ex x + —a) + Yd_ 



Making again the ansatz 

ip(x) 

and defining 



if)(x) = X^oipix) 



VP 



<p(xi) 



. it 2ira 



we arrive at the coupled set of equations 



d 



Ey + — \if 
v ay J 

E y ~ -J-] <P 



-\ip 
-\{p 



where we have set y = x m . Define the operators 



V2 



E\y + 



E\y 



1 d 

f \E\dy_ 

1 d 



(3.19) 



(3.20) 



(3.21) 



(3.22) 



(3.23) 
(3.24) 



Vz \ v 1 K ' xf^dy, 

These operators evidently satisfy the commutation relations of destruction 
and creation operators, respectively: 

[a, a'] = 1 



Substituting one equation into the other in (|3.22| ) we have, depending on the 
sign of E, 

2\E\aa [ <p = \ 2 tp 



E positive : 



2\E\a)aip = \ 2 <p 



(3.25) 
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and 

( 2\E\da(p = \ 2 (p 
E negative : < \ (3.26) 

[ 2|S|aoV = 

Now, 2|£ , |a^a is just the Hamiltonian of the harmonic oscillator with 
the zero-point energy ommitted. Correspondingly ip and (p are given by 
the harmonic oscillator eigenfunctions. Defining the ground state |0 > by 
a\Q >= we conclude that the eigenstates and corresponding eigenvalues are 
given by 



E positive: |tf(±) >= , \ n = ± J2n\E 




X n = ±y/2n\E\ 

where \n > are the eigenstates of the harmonic oscillator and are the 
corresponding energy-eigenvalues without the "zero-point energy" . Denoting 
by (p n (xi) the eigenfunctions of the harmonic oscillator, normalized with 
respect to the interval [—00, 00] and setting y = x m , we have for n > 1 the 
orthonormalized eigenfunctions of the Dirac operator ( 3-19Q , 



O) = ^ i(2m+1) **° 



y T¥n-l{Xr. 



n > 1 , (3.27) 



for positive E, and 




Wm,n( X ) = ^/^ e \ ,n>l, (S.2b) 



for negative E, each corresponding to the eigenvalues 



X n = ±\]2n\E\ , (3.29) 

respectively. Since the spectrum corresponds to the absence of the "zero- 
point energy" of the harmonic oscillator, we have an infinite set of orthonor- 
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malized zero modes labelled by m and chirality, of the form 

/ 



,(+)/ \ 1 i(2m+l)^x° 

0mW = ^e v >P 



V o 



(3.30) 



for positive E, and 

/ 



(x) = 7T 



(3.31) 





for negative E, each corresponding to the eigenvalue Ao = 0. This is in line 
with the Atiyah-Singer Index theorem in the infinite volume limit (see |J). 
Notice that in the case of the zero-modes, the superscript denotes " chirality" . 
The wave functions ( |3.30| ) and (|3.31| ) correspond to the ground state wave 



functions of the harmonic oscillator, localized at the positions X\ = (~"' +J '" 



^ with m G Z . This provides a physical interpretation of the degeneracy of 
the spectrum. In order to gain a further insight into the problem, we examine 
next the effective Lagrangian giving rise to this degeneracy, as defined in 
terms of the "local" ^-function. 

In order to simplify the discussion, we shall restrict ourselves in the fol- 
lowing to the case where E is positive. 

Effective Lagrangian density 

We begin by considering the local heat kernel. For the case in question 
it takes the form (we now take E > 0; we include the zero modes.) 

oo / oo 



h { S(t; x,y)= £ EE z- 2nEt ^U x )^Uv)e + fflWa&Xv)} 

m=— oo \n=l cr=± 

or explicitly 



OO 1 

hj${t;x,y) = £ L e ^n^-y°) (3.32) 



m=— oo 



En= e- 2nEt Vn(XmM(ym) ^ 

££ie-^»-i(s m K-i(j/m) ) 
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The diagonal matrix structure is again a consequence of the existence of 



a pair of eigenf unctions tp!^' corresponding to the eigenvalues ±V2nE, if 
We now observe that (note that the sum starts with n = 0) 



n=0 



-2nEt 



Ti 



(3.33) 



where the matrix element on the r.h.s. is the propagation kernel of the 
harmonic oscillator known to be given by (] 



< x\e- tHuo \y >- 



-Et 



7T s/l 



-4Et 



cxp 



E (x 2 + y 2 ){l + e~ AEt ) - Axye 



2Et 



1-e 



-AEt 



(3.34) 



Going to the limit of coincident points x = y, and taking the trace in matrix 
space, we arrive at 



trh%l(t;x,x) 



OO 

£ - [2coshEt < x m \e- tHHO \y m > 

m=—oo 

1 [e I 



= V e 

/3 V 7T Vtanh m= _oo 



-Ext, tanh St 



(3.35) 



Making use of the identity ( 2.34 ) we may thus write the heat kernel (|3.35| ) 
in the form 



tr/i (/3) (t 



E ( 



1 + 2 J2 (~ l T cos [m{E(3xi + 2vra)] e~ ™ft 



m=l 



2tt VtanhM, 

In order to compute the effective Lagrangian density we first need to subtract 
the zero-mode contribution: 



trti {/3 \t;x,x) 



tr/i (/3) (t;x,x) - - ^ (Po(Xm)<Po(%m) 
m=— oo 



where y^o^) is the zero-energy harmonic oscillator wave function: 

ip (x) - 



E\* _E2 

— I e 2 X 

7T 



Using fl2.3|) we have 
tr/i /(/3) (t;x,x) = — 

27T 



/o(t) + 2 £ (-l) m cos [m{Ef3x x + 2vra)] / m (t) 



(3.36) 

9 the Hamiltonian in our case is of the form H = p 2 + E 2 y 2 , and thus correponds to 
making the identifications m = i, u> = 2E in the conventional hamiltonian. 
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where 

fm(t) 



g 4tanh£t — 

tanh Et 



where the "prime" indicates the exclusion of zero-modes. From here we 
obtain for the effective Lagrangian 



^C (/3) (w) 



+ ^C(0;x,x) ln/i 2 
s =o 1 



where 



(^\s-x,x) = — / dtt s - l ti^\t-x,x) 
1 s JO 



and /i is an arbitrary scale paramenter reflecting the usual ambiguity asso- 
ciated with a change in scale of the dimensionful eigenvalues A n . We see 
that £ e ff(xi) is again a periodic function of x% with period AT = J^. The 
degeneracy of the A n spectrum with respect to the Matsubara index labelling 
the zero modes, reflects this fact. It is clear that correspondingly the energy 
density will exhibit the same periodicity. Notice that we have again obtained 
a clean separation of the /^-independent (sea) and /3-dependent (plasma) con- 
tributions. 

Effective action 

Restricting the dimension of our system to k "potential wells" (k zero 
modes) we have the following relation between E(3 and the length L of our 
system: 

L = k^. (3.37) 

Integrating C e jf over a space-time volume (3L we thus have for the effective 
action 

5e// = ^[C , (0) + C(0)ln / , 2 ] , (3.38) 

10 In the ^-function regularization the ambiguity in the calculation of the effective action 
is well known to be determined by C(0): In AetA — — C'(0) + £(0)ln/i 2 , where (x is an 
arbitrary scale parameter. 
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where 

C(s) = / dx I dx!C m (s;x,x) . (3.39) 
Jo 

and the "prime" now means differentiation with respect to s. From (|3.36| ) 
we have (the cosine term does not contribute to the integral) 

k r°° 

coo = L dtt s - i f (t) 

s-l e ~2nEt 



T(s) Jo 
2k 00 



00 1 9k 

2A: y t— -r = t^Cr(s) (3.40) 
^ (2n£) s (2£)* SKV ; v ; 



where (r{s) is the Riemann (^-function 



00 1 

Us) = Ezi 



n=l ™ 



Differentiating with respect to s, setting s = and using Cr(0) = — \, C'(0) 
' ln27r, as well as ( |3.37| ), we obtain for the effective action 

5 e// = -lnZ = ^ln(-^) (3.41) 

in agreement with the result obtained for the corresponding functional de- 
terminant on the torus in the presence of a finite number of zero modes || . 
From (|3.41| ) we have for the average thermal energy 

U=-^-lnZ = ^H^- 2 ). (3.42) 
op Air ti/i 2 

Since the energy is temperature independent there is no way of normaliz- 
ing it relative to the T = case. Notice that this temperature independence 
of the total thermal energy is a consequence of the spacial periodicity of C e ff, 
which in turn is the result of gauge invariance and periodicity in the time 
direction. 



4 Conclusion 



We have investigated the effect that uniform background electric and (in 
less detail) magnetic fields have on the distribution of thermal matter fields 
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in unbounded space. By calculating the thermal stress tensor Tj£>(x) and 
effective Lagrangian density we found that in the presence of a constant 
electric field the thermal plasma distribution of both scalar and fermion fields 
becomes periodic along the direction of E with period Ax = while the 
virtual sea remains uniform. On the other hand, a background uniform 
magnetic field does not lead to a spatial non-uniformity in either plasma or 
sea unless we compactify one of the spacial directions perpendicular to the 
magnetic field (say x 2 ) to a circle. In this case periodicity of the energy 
density distribution is again obtained along a direction x\ perpendicular to 
both 2 2 and to the magnetic field. In both cases (electric and magnetic field) 
this periodicity can be traced to the gauge invariance of the theory and the 
possibility of mapping Aq{x) and A 2 (x) into the intervals [0, and [0, ^], 
respectively, by a bonafide gauge transformation. This periodic structure 
was shown to reflect an infinite degeneracy of the eigenvalue spectrum of the 
spectral operator. 

It is interesting to note that the case of constant magnetic field in 3 spatial 
dimensions with one of the spacial dimensions in the plane orthogonal to the 
magnetic field compactified, i.e., the quantum Hall problem on a cylinder, 
is formally equivalent to the finite temperature case with constant electric 
field in 1+1 dimensions. This stems from the fact that at finite temperature 
we are working in 2 dimensional Euclidean space with temporal direction 
compactified. In the latter case the degeneracy of the Landau levels in the 
quantum Hall Effect corresponds to the number of zero modes of the Dirac 
operator as implied by the Atiyah-Singer Index theorem. 

Our exact results for a constant background electromagnetic field have 
a natural extension to arbitrary static background fields E(x) and B(x). 
One relevant equation to consult is eq. (|2.9|) . There we see that the space 
components of the vector potential are decoupled from the Matsubara index 
m. For an arbitrary static magnetic field B = V x A with A = (recall 
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that at finite temperature Aq = cannot be generally achieved by a gauge 
transformation) one simply makes the replacements }? mn — > (^p) 2 and 
<fmn ¥n in eq. ( |2.25| ) . The spatial mode equation then becomes 



-{V-iAf 



For arbitrary B(x) one thus has a situation much like the non-gauge theory 
of section 2.1, leading to some non-uniform distribution of both the sea and 
plasma components comparable to what physical boundaries would cause. 

The situation is quite different when A Q is non-zero. Then the static 
background field affects the sea and plasma quite differently. While plasma 
periodicity along the direction of E is strictly true only for constant E, one 
would expect for any electric field which is only weakly dependent on x 
a roughly periodic response from the plasma, and a nearly uniform distri- 
bution of the sea-component. Mathematically we have an infinite set of 
eigenvalue equations with a different Schrodinger like background poten- 
tial V m (x) = [Aq(x) — 27rm/j3} 2 for each Matsubara frequency. Note that 
m — > m + iV corresponds to performing an allowed gauge transformation 
with gauge function A = x q (2ttN/ f3), and that V m+ N(x) and V m (x) are con- 
nected by this gauge transformation. This situation differs fundamentally 
from the non-gauge case, where the potential V(x) does not bear the label 
m. Since Green functions, the energy momentum tensor, etc., are given in 
terms of equally weighted sums over all the individual problems labelled by 
m, they are explicitely gauge invariant. 

In general the diagonal heat kernel of a scalar or fermion quantum field 
at finite temperature T > in a static background is expected to factorize 
in the following way: 

c) = h(t; x, x) T =o [1 + fit; x; T)] 

where h(t; x, x)t=o is the temperature- zero heat kernel for the same back- 
ground, and f(t, x, T) is some function of the temperature T, the diffusion or 
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"proper" time t and the spatial position x. This function / vanishes exponen- 
tially as either T — > or t — > 0. The factorization above is motivated by the 
expectation that h^(t; x, x) separates quite generally for a static background 
into an UV divergent sea part, and an UV finite gas part 

^ (t, X, x) ^)sea ~t~ ^(^7 %)gas 

where 

x, x) sea = /i(t; x, x) T =o • 

Defining f(t;x;T) by 

^) 9 a S = fit; x] T)h(t; x, x) T =o 

we arrive at the factorization above. Known properties of h(t; x) gas lead to 
the stated properties of /. 

It is a matter of some interest to study the function /(£, x; T). Let us list 
the explicit examples computed in this paper. 

Nongauge scalar theory: 



i + /= E 



e 



n 2 f3 2 /4t 



Gauged scalar theory with A = Ex\ + 2ira/j3: 



i + / = E 



e -n 2 f3 2 E/4tanhEt e -inl3A 



Schwinger model with A = Ex\ + 2na/f3: 



i + / = E 



_ \ n — n/3 2 E/4 tanh Et — in(3A 



n=— oo 



From the gauge theory examples we see that f(t, x; T) in general depends on 
when a gauge potential background is involved. 
In ref. fl2| it is argued that / has the simple form 



i + /= E 



for a gauged spinor theory in a general (even time- dependent) background 
gauge potential A^. This claim is incorrect: the result for the Schwinger 
model provides a simple, explicit counterexample. 
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